Double jumps and transition rates for two dipole-interacting atoms 
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Cooperative effects in the fluorescence of two dipole-interacting atoms, with macroscopic quantum 
jumps (hght and dark periods), are investigated. The transition rates between different intensity 
periods are calculated in closed form and are used to determine the rates of double jumps between 
periods of double intensity and dark periods, the mean duration of the three intensity periods and 
the mean rate of their occurrence. We predict, to our knowledge for the flrst time, cooperative effects 
for double jumps, for atomic distances from one and to ten wave lengths of the strong transition. 
The double jump rate, as a function of the atomic distance, can show oscillations of up to 30% 
at distances of about a wave length, and oscillations are still noticeable at a distance of ten wave 
lengths. The cooperative effects of the quantities and their characteristic behavior turn out to be 
strongly dependent on the laser detuning. 

PACS numbers: 42.50.Ar, 42.50.Fx 



I. INTRODUCTION 

The dipole-dipole interaction between two atoms can 
be understood through the exchange of virtual photons 
and depends on the transition dipole moment of the lev- 
els involved. It can be characterized by complex coupling 
constants, or by their real and imaginary parts, where the 
former affect decay constants and the latter lead to level 
shifts [H . Cooperative effects in the radiative behavior of 
atoms which may arise from their mutual dipole-dipole 
interaction have attracted considerable interest in the lit- 
erature Two of the present authors have in- 
vestigated in detail the transition from anti-bunching to 
bunching with decreasing atomic distance for two dipole- 
dipole interacting two-level atoms. 

The striking phenomenon of macroscopic quantum 
jumps (electron shelving or macroscopic dark and light 
periods) can occur for a multi-level system where the 
electron is essentially shelved for seconds or even minutes 
in a metastable state without photon emissions 
For two such systems the fluorescence behavior would, 
without cooperative effects, be just the sum of the sepa- 
rate photon emissions, with dark periods of both atoms, 
light periods of a single atom and of two atoms. In Ref. 
pQI the fluorescence intensity of three such ions in a Paul 
trap was measured and a large fraction of double or triple 
jumps was reported, i.e. jumps by two or three intensity 
steps within the short resolution time. This fraction was 
orders of magnitudes larger than that expected for inde- 
pendent ions. A quantitative explanation of such a large 
cooperative effect for distances of the order of ten wave 
lengths of the strong transition has been found to be dif- 
ficult H, H, ^, 0. Other experiments at larger 
distances and with different ions showed no cooperative 
effects H, |§. 

Quite recently, two of the present authors inves- 
tigated for two such systems cooperative effects in the 



mean duration, Tqi Ti, and T2, of the dark, single- 
intensity, and double-intensity periods, respectively. This 
was done by simulations for two atoms in a configu- 
ration. The mean duration of the single- and double- 
intensity periods depended sensitively on the dipole- 
dipole interaction and thus on the atomic distance r. 
They exhibited noticeable oscillations which decreased in 
amplitude when r increased. These oscillations seemed 
to continue up to a distance of well over five wave lengths 
of the strong transition and they were opposite in phase 
with those of ReC3(r), where C3 is the complex dipole- 
dipole coupling constant associated with the strong tran- 
sitions. 

In this paper we present an analytic approach to study 
cooperative efl:ects for atoms with macroscopic quantum 
jumps. This is explained for two atoms, but is easily gen- 
eralized. The approach is based on an explicit calculation 
of transition rates between the various intensity periods. 
From the transition rates all interesting statistical quan- 
tities can be determined, such as double jump rates and 
mean duration of different periods. 

We predict, to our knowledge for the first time, coop- 
erative effects in the double jumps of two dipole-dipole 
interacting atoms. These results are for atoms in the V 
configuration (see Fig. |l|) and are verified by simulations. 
As a function of the atomic distance, the double jump 
rates show marked oscillations, with a maximal differ- 
ence of up to 30%, decreasing as I/r. Most surprising is 
a change in the oscillatory behavior of the double jump 
rate from in phase with Re C3 (?') to opposite in phase 
when the detuning of the laser driving the weak atomic 
transition is increased. For the mean durations Ti and T2 
there can be a change in behavior from opposite in phase 
to in phase with ReC3(r). Moreover, for a particular 
value of the detuning, which depends on the other pa- 
rameters, the double jump rate becomes constant in r and 
the cooperative effects disappear. This is true also for the 
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mean period durations and for their mean occurrences, 
with different values of the detuning, though. Typically, 
for nonzero detuning the oscillation amplitudes do not 
exceed those found for zero detuning. 

The experiments of Ref. exhibited extremely large 
cooperative effects, in fact up to three orders of magni- 
tude. Since this was for a different atomic level configu- 
ration and for three ions in a trap our results do not apply 
directly. In principle, however, our analytic approach can 
be carried over to the experimental situation of Ref. , 
although the calculations become algebraically more in- 
volved and have not been carried out so far. 

The plan of the paper is as follows. In Section II the 
fluorescence with its three different intensity periods is 
treated as a three-step telegraph process and the Bloch 
equations are used to derive the transition rates between 
the periods. In Sections III and IV expressions for the 



double jump rate and the mean duration of the three 
types of intensity periods are obtained by means of these 
transition rates. The results are compared with simu- 
lations in which photon intensities are obtained by av- 
eraging photon numbers over a small time window. It 
turns out that this data-smoothing procedure can affect 
the results, and we show how this can be corrected for 
quantitatively. A similar effect can also occur when pho- 
ton detectors measure the intensity of light by averaging 
over a small time window. In the last section the re- 
sults are discussed. It is suggested that the mean rate of 
double-intensity periods is an experimentally more eas- 
ily accessible candidate for exhibiting cooperative effects 
arising from the dipole-dipole interaction. 



assume the relations 

r!2«r!3, n2^nl/A3, A2«0, A3-O, (1) 

A2 arbitrary. The Dicke states are defined as 

\g) = |1)|1), |e2) = |2)|2), leg) = |3)|3) 
\s,k) - {\j)\k) + \j)\k)}/V2 
= {b->|fc)-b-)|A:>}/V2 

They are symmetric and antisymmetric, respectively, un- 
der permutation of the two atoms. The Dicke states and 
the possible transitions are displayed in Fig. ||. Solid 




FIG. 2: Dicke states. Simple arrows denote decays. Solid 
and dashed double arrows denote strong and weak driving, 
respectively. 



II. TRANSITION RATES 
A. Prerequisites 

We consider two atoms, at a fixed distance r, each a 
V configuration as shown in Fig. nl We assume the 
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strong laser field 




weak field, fl-z 



FIG. 1: V system with metastable level 2 and Einstein coef- 
ficient A3 for level 3. Q2 and ^3 are the Rabi frequencies of 
the two lasers driving the weak 1-2 transition and the strong 
1-3 transition, respectively. 

laser radiation normal to this line, and for the Einstein 
coefficients, the Rabi frequencies and the detuning we 



single and double arrows indicate decay and strong driv- 
ing by laser 3, respectively, while dashed double arrows 
indicate the weak driving by laser 2. For — 0, i.e. 
with the dashed arrows absent, the states decompose into 
three non-connected subsets, namely |e2), the four states 
of the inner ring and the four states of the outer ring in 
Fig. H, and the subspaces spanned by these states will be 
denoted by dark, inner, and outer subspace, respectively. 
As in Ref. [|7| they will be associated in the following 
with the fluorescence periods of intensity 0, 1, and 2: 

dark state : |e2) (2) 
inner states (intensity 1) : |si2), |s23), |ai2), |a23) (3) 
outer states (intensity 2) : \g), jsis), jea), jaia) (4) 

The weak laser will lead to slow transitions between the 
subspaces. 

The Bloch equations can be written, with the condi- 
tional Hamiltonian iJcond and the reset operation TZ of 
Appendix in the compact form |5^ 

P=-J- [^^condP - pHl^J + n{p) . (5) 

The operator iJcond is of the form 

i^cond = -f^cond + -f^cond(^2) (6) 
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where the operator -ff^ond depends on fl^ and on the 
dipole-dipole couphng constant C3(r), while H^ond 
ear in ^2 and does not depend on C3 and fl^ . The super- 
operator TZ depends on C3. 



B. Intensity periods and subspaces 

For a single atom as in Fig. with macroscopic hght 
and dark periods, the stochastic sequence of individual 
photon emissions can be directly analyzed by the quan- 
tum jump approach g |9|, |o[ |l|, |2[ |5|, using 
the existence of different time scales. To high precision 
it yields a telegraph process and the transition rates be- 
tween the periods p^ , p6| . A more heuristic approach 
assumes that during a light period the density matrix 
of the atom lies in the subspace spanned by |1) and |3) 
and that during a dark period the state is given by |2) 
psf . One can then use the Bloch equations to calculate 
the build-up, during a time At, of a population outside 
the respective subspace and obtains from this the prob- 
ability of leaving the subspace. This probability is then 
interpreted as the transition probability from one period 
to the other. The results agree with those of the more 
microscopic quantum jump approach ||3^, |4^, [5^ . 

This idea will be used here for two dipole-interacting 
V systems. We associate each of the three types of flu- 
orescence periods with one of the subspaces spanned by 
the states in Eq. (H) - (^) and model transitions be- 
tween periods as transitions between the corresponding 
subspaces. Without dipole interaction this is the same 
assumption as for a single atom, and with the interac- 
tion it has been tested numerically in Ref. |Q to hold 
as long as the atomic separation is larger than a third 
wavelength of the strong transition. 

Thus, at a particular time to, the density matrix p{to) 
of the two atoms is assumed to lie in one of the subspaces. 
Then, during a short time At, satisfying 



n-\A-^ < At < n-\ 



(7) 



the system will go over to a density matrix p(tQ + At) 
which contains small populations in the other subspaces, 
due to the driving by ^ 0. The time derivatives of 
these populations at to + At give the transition rates to 
these subspaces because, as shown in Appendix B, they 
are independent of the particular choice of At and of 
the particular density matrix p{to), as long as Eq. (0) 
is fulfilled. These rates can be interpreted as transition 
rates between corresponding intensity periods, just as in 
the one-atom case. 

A straightforward calculation using Eq. (|) yields the 
exact relations 

^ ^ (outer |p|outer) 



flolm 



2(si2|p|5) + (S23|p|si3) + (a23|p|ai3)|8) 



^(e2|p|e2) = V2fi2lm(si2|p|e2) 



(9) 



dt ^ 

inner 



[inner I/) I inner) 



-^{(e2|p|e2) + (outer|p|outer)](10) 



where jouter) stands for \g), jsia), jes), |ai3) and jinner) 
for |si2), |s23), |ai2), |a23)- Thus one has to calculate the 
coherences on the right-hand side at time to + At to first 
order in ^2 , with the appropriate initial condition at time 
to, to obtain the transition rate to second order in f22- 

If p{to) lies in one of the above subspaces then by time 
to + At the system has reached a quasi-stationary state 
satisfying 



p{to + At) = to first order in il2, 



(11) 



as shown in Appendix B. This is also true for a single 
atom and is the decisive equation. 

To obtain from this the coherences to first order in 02 
we write 

p(to + At) = p" + pi + • • • 

where p'^ is of order fif. Putting p = in Eq. (j^) and 
inserting the expansion for p one obtains in zeroth order 



= 



i/,VdP°-P°<ld +np') (12) 



and in first order in 



-'^condP - P ^cond 



Wl ^0 ^Orrlt 
■'^condP - P ^con. 



-7^(pl)(13) 



Thus p° is an equilibrium state for = 0, taken to lie 
in the appropriate subspace. For the dark state and the 
subspace spanned by the inner states one has 



n° = n° 
P — Pdark 



Ie2)(e2| 



(14) 



^°^P?n„e..-^{pi^^®|2)(2| + |2)(2|0p(f)} 



1 + ^1"^ 

4^rq:^{ki2>(si2| + |fli2)(ai2|} 

1 n'^ 

^ {|S23)(S23| + |a23)(a23|} 



^^3^3 



2 ~\~ 2i^^ 



4:Al + 2nl 

{|S12)(S23| - |ai2)(a23|} + H.C. (15) 



by symmetry, independently of C3, where pit'^^ are the 
steady states of the individual atoms in the 1-3 subspace 
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(for ^2 = and C3 = 0). For the subspace spanned by 
the outer states one calculates 



0° = 0° 

P — Pouter 



{iAl + nir + Al\Cs\' + 2AlReC,}\g){g\ 

+ {iV2A3 rig + nt + A3 C3)|5)(si3l + H.c.} 
-{Aal)^ (A3 + C3)|ff)(e3|+H.c.} 

(2^2 + r!2)|si3)(si3| + r!t {|e3)(e3| + |ai3)(ai3|} 



we write {\xi)} = {|si2), IS23), |ai2), |a23)} (inner states) 
and {\yj)} = {\g), |si3), |e3), |ai3)} (outer states) for the 
corresponding bases and decompose 



P 



P^j I + PI3 \yo)i^^ + otlier terms (17) 



-{^V2A3 r!3|si3)(e3|+H.c.} 



(16) 



Inserting this into Eq. (^3|) and taking matrix elements 
with (xio I on the left and \yja) on the right gives 



One checks that for C3 = this becomes Poster ^ Pi' 
Pss , the expression for two independent atoms. 

We will denote the transition rates between the sub- 
spaces by pij. Here i,j = 0, 1, 2 refer to the dark, inner 
and outer subspace, respectively, (and thus to the cor- 
responding intensities). The pij will be determined to 
second order in ^2- As expected, po2 and P20 will turn 
out to be zero. 



C. Calculation of pi2 

We start from = Pilmcr ^1- (111) initial state. 
For the the transition rate pi2 to the outer subspace one 
needs, in view of Eq. (||), three coherences of be- 
tween the inner and outer subspace. To obtain these 



J 



= ^{Xia\P°nncvH^ond\y]o) - ^"^Pl^oi^^olHcc 



cond!'''*/ 



j 

+ ^(^3 + Re C3)pI^ {x,„ \R+\x,) {yj \Rl\yj„) 

+ ^(^3 -ReC3)p^,(x,Ji?_|x,)(y,|i?L|y,„) .(18) 

This is a system of 16 linear equations for the 16 coher- 
ences p]j, of which only three are needed in Eq. (||). 
Due to the symmetry of ffcond and ii+ and antisymme- 
try of i?_ under the interchange of the two atoms, the 
system decouples. Taking for \xig) and \yjo) either both 
symmetric or both antisymmetric states and putting the 
eight coherences into the column vector 



P — (Psi2g' /'si2Sl3' Psi2e3l Ps23S' Ps-23S13^ Ps23e3'' Pai2ai3^ Pa23ai3^ 



(19) 



one obtains the equation 

(A - iA2l)/5 = ai 



(20) 



where 



-i^s/V^ 
in3/2 














-in3/V2 

-in^/V^ {A3 + C^)/2 -1^3/ V2 



ai 



A3 



in3/2 



iD,2^3 



10.3/2 




A3/2 
-1^3/^2 






-(A3 + ReC3)/V2 

in3/2 


-1^3/^2 
(A3+C3V2) 
1^3/^/2 





.+ReC3)/V2 
10.3/2 


-^3/^/2 
3A3/2 



—I 



_(A3-ReC3)/\/2 {A3-C. 

n "■'^ 



A{Al + 201) 



V2^ 



AI 



O3 



, i ^3, 0, -i\/2 A3, 1^3, 0, -i A3, 1^3 



:)/2 



-103/2 



(21) 

(ReC3-A3)/V2- 






~i03/2 
(A3-C3V2) . 

(22) 



r 



Inverting the 8x8 matrix A — iA2l by Maple yields p illuminating. Inserting the required coherences into Eq. 
and the coherences. The result is complicated and not 
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one obtains, to first order in Re C3 and Im C3 and to 
second order in 0,2, 



P12 



n. 



+ReC3{r) 



2Alnl{nj-4:AlAl-16A 



{Ai + 2ni) {nj - sA^nj + aa^ai - 



-16A4)2 



}23) 



Note that only Re C3 appears and that the terms linear 
in Im C3 have canceled. 



D. Calculation of pio 



To determine pio we use Eq. and start again from 
p° = /Conner initial condition in Eq. (p^), but now have 



to determine 
choosing 



vSl2|P 162; 
= |S12), 



Replacing {\i 
\s23) in Eq. m 



} by |e2) and 
I), one obtains 



two inhomogeneous linear equations for (si2|yO |e2) and 
(s23|p^|e2). These equations do not depend on C3, since 
R± and vanish on |e2) and since C3 does not appear 
in the part of -ffcond acting on the inner states. Therefore 
(si2|p^|e2) and (s23|/9^|e2) are independent of C3. By a 
simple calculation one obtains (si2|p^|e2) and inserting 
this into Eq. (0) yields, to second order in 0,2, 



to start from p° = plj^rk initial condition in Eq. jl^). 
Therefore one obtains the same equations for (si2|p^|e2) 
and (s23|p^|e2) as before, except for the inhomogeneous 
part. One has independence of C3 and easily solves for 
(si2|p^|e2). 

For the remaining coherences needed in Eq. ([l0|), i.e. 
those in Eq. (^), one obtains the same form as in Eq. 
(pO|), with the same A, but now with ai = since the 
term containing p'^ vanishes. Therefore these coherences 
vanish here and hence po2 = 0. Physically this means 
that in our formulation of the problem there are no direct 
transitions from a dark period to a period of intensity 2. 



From Eq. (hOl) one now obtains, to second order in fl2 



Poi = 2nl 



A^nl 



AAlAl 



(25) 



This is independent of C3 and is the same as for two 
independent atoms, namely twice the transition rate for 
a single atom from a dark to a light period. 

F. Calculation of p2i 



PlO 



A^nliAl 



-4Ai) 



{Ai + 2ni)mi-iAiy + 4AiAi] 



(24) 



This is independent of C3 and is the same as for two in- 
dependent atoms, namely the transition rate for a single 
atom from a light to a dark period [E9[ . 



E. Calculation of poi 

To determine poi we use Eq. (p^. One also needs 
(si2|/0^|e2), as seen from Eq. W), but in this case one has 



The transition rate p2i is obtained from Eq. ( |10[ ) and 
the required coherences are again those appearing in Eq. 
(^ and (p|), now with = Poutor ^ initial condition. 
For (si2|p |e2) and (s23|p^|e2) one obtains the same two 
equations as before, except for the inhomogeneous part 
which now vanishes. Hence these two coherences vanish 
now and as a consequence P20 = 0. Physically this means 
that in our formulation there are no direct transitions 
from a period of intensity 2 to a dark period. For the 
coherences in Eq. (^9|) one has the same equation as Eq. 
(P0|), with the same matrix A but with ai replaced by 



a2 = -1^2 /{ V2 (4 1^3^ + 4 n^^As^ + A^^Re C3' + 2 A3'Rc C3 + A3'lm C3' + A3*) } 



(26) 



+ 2 na^A^' + As'Re Cy + 2 A3^Re C3 + A3'^Im C3' 
ifl^ V2A3 (A3' + A3 Re C3 + ilm C3 A3 + ^3^) 

-^83' (A3 + Re C3 + ilm C3) A3 
in3 A3 {~n3^ - A32 - A3 Re C3 + lira C3 A3) 
n3^n3^ + 2A3^) /V2 

iri3'A3 



f^3V\/2 



A3 
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Inserting the resulting coherences into Eq. ( jioj ) gives, to 
first order in Re C3 and Im C3 and to second order in j 
I 



^2; 2^3171(^1 + 4A 



4A2l]2(^4^4 ^ 4^2^6 _ i2A2A2r!^ - UAjl^l - 4A6A2 - 32A4A4 ^ 64A|1]| + leAjlll) ^ 

^^"■^ (A2+2f]2)3(04_8^2^2_^4^2^2^1g^4)2 / 



(28) 



where again the terms containing Im C3 have canceled. 

G. Discussion. 

If one computes the coherences in Eqs. (||) and ( |lO| ) 
to second order in C3 one obtains pi2 and p2i to second 
order in C3. The resulting expressions are not enlight- 
ening and therefore not given here, but they do depend 
on (ImC3)2. Fig. || shows how small the second-order 
dipole-dipole contribution to p2i is for the parameters of 
the simulations and for distances larger than half a wave 
length. For smaller distances the results are probably not 



P12 and P21 vary with the atomic distance in phase with 
ReC3. For A2^0, however, the coefficients of ReC3 in 
Eqs. (H^) or (g^) can become zero or negative. In the 
first case pi2 or p2i become constant in r, while in the 
second case they vary opposite in phase to Re C3 . 

It will be shown in the next sections that this depen- 
dence of pi2 and P21 on the detuning of the weak laser 
entails a corresponding behavior of the double jump rate 
and an opposite behavior of the mean durations Ti and 
T2. This opposite behavior of Ti and T2 is easy to un- 
derstand since they are related to the inverse of the tran- 
sition rates. 



P21 to 1st order in C3 ■ 
P21 to 2nd order in C'3 
P\2 to 1st order in C3 " 




30 
tail' 



FIG. 3: Transition probabilities p2i to first and second order 
in C3 and p2i to first order, for 173 = 0.5^3, Q.2=G.QlAz, zero 
detuning. Tlie contribution to P21 arising from the second 
order in C3 is small. 

applicable anyway, as discussed in Ref. p7| . 
For A2 = the rates pi2 and p2i simplify to 



Pi2 = r!^{^+RcC3(r) 



2AI 



} (29) 



P21 



2Al 



^3(^3 



ReC3(r) 3^ 3 3' 



f]2(Ai + 2f]2)3 



}(30) 



and one sees that the coefficients of the ReC3 term in 
Eqs. (|2|) and ( pO[ ) are positive. For A2 = 0, therefore. 



III. DOUBLE JUMPS: COMPARISON OF 
SIMULATIONS WITH THEORY 

A double jump is defined as a transition from a double- 
intensity period to dark period, or vice versa, within a 
prescribed time interval ATdj. Now, to distinguish dif- 
ferent periods in experiments and in simulations one has 
to use an average photon intensity, obtained e.g. by 
means of averaging over a time window. This window 
has to be large enough to contain enough emissions, but 
must not be too large in order not to overlook too many 
short periods. Our simulations employ a procedure sim- 
ilar to that in Ref. |4^ and use a moving window of 
fixed width, denoted by ATw The time interval ATdj 
should be larger than ATw. 

We consider the fluorescence periods as a telegraph 
process with three steps and use the pij of the last section 
as transition rates. At first the influence of the averaging 
window r„ will be neglected. 

The rate of downward double jumps is obtained as fol- 
lows. For i = 0, 1, 2, let be the mean number of 
periods of intensity i per unit time. For a long path of 
length T the total number of periods of intensity i is then 
Ni(T) = UiT. At the end of each period of intensity 2 
there begins a period of intensity 1, and the probability 
for this period of intensity 1 to be shorter than ATdj is 
given by 

1 - exp{-(pio +pi2)ArDj}. 

At the end of a period of intensity 1 the branching ratio 
for a transition to a period of intensity ispio/(pio+Pi2)- 
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Thus during time T the total number of such downward 
double jumps, denoted by N^j{T), is 



N^",{T)^N2iT) 



Pw 



(Pl0+Pl2 



■|l-exp{-(pio+pi2)ArDj}| 



and therefore the rate, n^j, of downward double jumps 
within ATdj is 



"^DJ="2 _^ {l~ exp{-(pio +Pi2)ATdj}|. 



(P10+P12) 



(31) 



In a similar way one finds that the rate, n'^j, of upward 
double jumps within AT^j is 



02 P12 



(Pl0+?'12^ 



exp{-(pio +pi2)ATdj}|. 



(32) 



It remains to determine uq and 712- Since a period of 
intensity 1 ends with a transition to a period of either 
intensity or intensity 2 one has, with the respective 
branching ratios, 



Pw 

no = ; ni 

P1Q+P12 



P12 

712 = "1- 



Pio +P12 



(33) 



(34) 



If one denotes by the mean durations of a period of 
intensity i, one has 



(35) 



Moreover, one has 

To = l/poi, ^1 = l/(pio +P12), 72 = 1/^21 (36) 
and this then gives 

P01P21 



no 



n2 



P01P21 +P21P10 +P01P12 
P01P21 



-Pio 



-P21- 



(37) 



(38) 



P01P21 +P21P10 +P01P12 

From this, together with Eqs. (|l|) and (||), one sees im- 
mediately that the rates of upward and downward double 
jumps are equal. 



„02 _ 20 



(39) 



This fact was also observed in the simulations. The com- 
bined number of double jumps therefore equals 



P01P10P12P21 



(P01P21 + P21PW + PoiPi2){poi + P12) 
X {1 - exp{-(pio +pi2)ATdj}}.(40) 



For ATdj <C Ti and by expanding the exponential, this 
gives for the combined double jump rate, without correc- 
tion for the averaging window. 



P01P10P12P21 



P01P21 + P21P10 + PQ1P12 



-ATdj 



(41) 



Fig. ^ shows a comparison of this result with data from 
the simulations. Except for atomic distances less than 




FIG. 4: Double jump rates. Simulation -I- + -I- , theory — 
— uncorrected for averaging window (Qa —0.5 A3, Q2 = 
0.01 A3, zero detuning). 

about three quarters of the wave length of the strong 
transition the agreement appears as quite reasonable, 
and the disagreement for small distances is not unex- 
pected since there the intensities start to decrease and 
a description by a telegraph process may be no longer a 
good approximation, as pointed out in Ref. But 
one observes that the theoretical result is systematically 
above the simulated curve. This seeming disagreement, 
however, is easily explained and can be taken care of as 
follows. 



A. Corrections for averaging window 

We recall that the simulated data were obtained by av- 
eraging the numerical photon emission times with a mov- 
ing window of length AT^. Then, roughly, periods which 
are shorter than about two thirds of the window length 
are overlooked, and therefore the number of recorded (or 
observed) periods of type 2, which enters Eq. (|3l|), is 
smaller than that given by Eq. ( |38| ) . The recorded or ob- 
served number is denoted by 712, cor- It is approximately 
given by 



«2.. 



7i2exp{-p2igATw}, 



(42) 



and this expression should be inserted into Eq. (^1|) for 
7X2. In this way one obtains the corrected theoretical 
curve m Fig. |. The curve changes very little if instead 
of two thirds one takes 60% or 70% of AT^. It is seen that 
the agreement with the simulated data is much improved 
for distances greater than three quarters of a wave length 
of the strong transition. 
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FIG. 5: As in Fig ^ , but theory corrected for averaging win- 
dow. 



It still appears, however, that the oscillation ampli- 
tudes of the theoretical curve are somewhat larger than 
those of the simulated curve. This is again understand- 
able as an effect of the averaging procedure. In the sim- 
ulations it was noticed numerically that the r depen- 
dence of the double jump rate depended somewhat on 
the length of the averaging window and distinct fea- 
tures tended to be somewhat washed out for larger AT^, 
in particular the oscillation amplitudes of the simulated 
data decreased with the length of the averaging window. 
A larger AT^ gave a smoother intensity curve, but made 
the determination of the transition times between dif- 
ferent periods more difficult, while a shorter averaging 
window introduced more noise. We found the use of 
ATw = 114 ^3^^ to be a good compromise. If it were 
possible to choose smaller averaging window the ampli- 
tudes should increase, as predicted by the theory. 



B. Detuning 

One can explicitly insert the expressions for pij of the 
last section into Eq. (pl|), but the result becomes un- 
wieldy. One can show that in an expansion of Eq. (^ ) 
with respect to Re C3 to first order the coefficient of Re C3 
is positive for zero detuning. This implies that the dou- 
ble jump rate is in phase with ReC3(r) for the atomic 
distances under consideration and for zero detuning. For 
increasing detuning the double jump rate can become 
constant in r and then change its oscillatory behavior to 
that of — ReCa. An example for the latter is shown in 
Fig. g 



IV. DURATION OF FLUORESCENCE 
PERIODS: EFFECT OF AVERAGING WINDOW 

The mean durations, Tq, Ti, and T2, of the three peri- 
ods were investigated for cooperative effects in Ref. Q 
by simulations with averaging windows at discrete times. 
Here we have performed similar simulations with a mov- 
ing averaging window. It turns out that both the present 
and the previous simulation for Ti are about 15% higher 
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FIG. 6: Changed oscillatory behavior of the double jump rate 
for increased detuning of the weak driving (uncorrected for 
averaging window). 



than those predicted by Eq. (|3q), using the expressions 
for pij of Section II and without correcting for the use of 
the averaging window due to which short periods are not 
recorded. We will now show how this can be taken into 
account in the theory. 



As in Section III we consider a three-step telegraph 



process with periods of type 0, 1, and 2, whose mean du- 
rations are denoted by Tq, Ti, and T2, respectively. We 
assume that periods of length At or less are not recorded. 
Fig. 1^ shows periods of type 1 which are interrupted by a 
short period of type and 2, respectively. If the respec- 



FIG. 7: If periods of length At or less are overlooked then 
the distribution of the periods is changed. 

five short periods are not recorded, then the two periods 
of type 1 in the left part of the figure are recorded as a 
single longer period, and similarly for the right part of the 
figure. This leads to an apparent decrease of shorter pe- 
riods of type 1 and to a corresponding increase of longer 
periods. 

To make this quantitative we put = 1 /Tt and denote 
the number per unit time of periods of type i, whose 
duration is less than At, by nf i.e. 



i|l -exp{-AiAr}}. 



(43) 



Per unit time, one has n^'^ occurrences of the situation in 
the left part of Fig. |^ and n^"^ occurrences of that in the 
right part. The probability for one of the periods of type 
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1 in the left or right part of Fig. ^ to have a length lying 
in the time interval (ii,ti + dti) is 2Ai exp{— Aiiijdii, 
where the factor of 2 comes from the two possible situa- 
tions. Therefore, the recorded number, per unit time, of 
periods of type 1 with duration in (ii, ti +dti) is changed 
(decreased) by 



2(n^^ + nf^)Xi exp{-Aiii}dti 



(44) 



Similarly, the apparent increase of the number, per unit 
time, of periods of type 1 with duration in (ti,ti + dti) 
is, by Fig. 0, 

(n^'+ntl J J dt[dt'{ 

ti<t[+t['<ti+dti 

X Ai exp{— Ai^ijAi exp{— Alt"} 
= (no'' + "^^)A?ti exp{-Aiti}(iti(45) 

Denoting by h'ircciti)dti the actually recorded number, 
per unit time, of periods of type 1 with duration in 
(ti,ii + dti) one obtains from the two previous expres- 
sions 

iyircc{ti)dti = niAi exp{-Aiti}dti 

+ {n^^ + n^^){Xlti ~ 2Ai) exp{-Aiti}dti(46) 

The average duration of the recorded periods of type 1 
will be denoted by Ti^cor, and it is given by 



Tl,cor — / dti ti {ti)/ / dtivi,,,{ti). (47) 

JAt I J At 

Using Eq. ( p6| ) for i^ircc(ii) one obtains, after an elemen- 
tary calculation and for Ar satisfying Ar/Ti <C 1, 



PW +Pl2 



^^^^^mPio+Pi2P2i (48) 



(PlO +Pl2)^ 



The first term is the ideal theoretical value, Ti, and the 
remainder is the correction due to non-recorded short 
periods. In a similar way one obtains 



Poi L Poi 



72,cor = h At <^ IH 

P2l L P21 



(49) 



(50) 



where again the respective first terms are the ideal values. 
To and T2 . 

To compare this with simulated data, obtained with a 
moving averaging window of length AT„ = 247 we 
have taken At = |ATw, as in the previous section, and 
have plotted the results together with the simulated data 
in Fig. ^ The agreement is very good. Quite generally, 
for zero detuning the oscillations of Ti and T2 are opposite 
in phase to those of Re G^{r), as already noted at the end 
of Section II. As in the case of the double jump rate, Ti 
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FIG. 8: Mean duration of fluorescence periods. Simulation: 

T2 + + +, Ti XXX, To ***. Theory: T2 , Ti — 

— , To .... corrected for averaging window (Sla=0.5 A^, 
r22 = 0.01 A3, zero detuning). 



and T2 can become constant in r for particular values of 
the detuning (different for Ti and T2 ), and then change 
to a behavior in phase with ReC3(r). 

The above approach of taking the averaging window 
into account works for the following reason. For a sin- 
gle atom with macroscopic dark periods it is known that 
the emission of photons is describable, to high accuracy, 
by an underlying two-step telegraph process. For two 
independent atoms with macroscopic dark periods the 
emissions are therefore described by an underlying three- 
step telegraph process. For two atoms interacting by a 
weak dipole-dipole interaction the actual emission pro- 
cess of photons should therefore still have, at least ap- 
proximately, an underlying three-step telegraph process. 
What we have done above is replacing the actual emission 
process by this underlying three-step telegraph process 
and then incorporating the averaging window by taking 
into account the influence of the overlooked short periods 
on the statistics. 



V. DISCUSSION OF RESULTS 

We have investigated cooperative effects in the fluo- 
rescence of two dipole-dipole interacting atoms vti a, V 
configuration. One of the excited states of the V configu- 
ration is assumed to be metastable, i.e. with a weak tran- 
sition to the ground state. When driven by two lasers, a 
single such configuration exhibits macroscopic dark pe- 
riods and periods of fixed intensity, like a two-step tele- 
graph process. A system of two such atoms exhibits three 
fluorescence types, i.e. dark periods and periods of single 
and double intensity, like a three-step telegraph process. 
For large atomic distances, when the dipole-dipole inter- 
action is negligible, the total fluorescence just consists of 
the sum of the individual atomic contributions. We have 
shown that for smaller atomic distances the fluorescence 
modified by the dipole-dipole interaction which depends 
on the atomic distance r. In particular we have, to our 
knowledge for the first time, explicitly demonstrated co- 
operative effects in the rate of double jumps from a period 
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of double intensity to a dark period or vice versa, both 
analytically and by simulations. 

By means of an analytical theory we have obtained the 
r-dependent transition rates , pij, between the three in- 
tensity periods. These were then used to calculate the 
rate of double jumps and in the mean period durations 
To, Ti, and T2. When comparing with the simulations it 
turned out that one had to take into account the averag- 
ing window used for obtaining an intensity curve from the 
individual photon emissions. With this the agreement be- 
tween simulation and analytic theory became excellent. 

For zero laser detuning, for which the simulations were 
performed, the double jump rates are in phase with and 
Ti and T2 opposite in phase to Re C^ir). The theoretical 
expressions, however, allow general detuning, A2, of the 
laser which drives the weak transition. It has been shown 
that for a particular A2 , which depends on the other pa- 
rameters, the double jump rate becomes constant and, 
for larger A2, varies opposite in phase to KeCslr). A 
similar change of characteristic behavior also occurs for 
Ti and T2, for different values of A2 though. The ampli- 
tude of the oscillations with the atomic distance remain 
in the same region of magnitude as for zero detuning. As 
pointed out in Rcf . Q , a dependence of the oscillations 
on Re C3 (r) is not unexpected since Re C3 (r) affects the 
decay rates of the excited Dicke states of the combined 
system. But an intuitive argument why the above change 
of behavior occurs for increased detuning is at present not 
apparent. 

that there is another 



We have pointed out in Section II] 



statistical property of the fluorescence which can serve as 
an indicator of the influence of the dipole-dipole interac- 
tion and which is probably not too difficult to determine 
experimentally. This quantity is the rate with which flu- 
orescence periods of definite type occur, in particular the 
rate of periods with double intensity. Our theoretical re- 
sults show that this rate behaves similar to the double 
jump rate, as regards the variation with the atomic dis- 
tance, and an example is shown in Fig. 0. This quantity 




rate or the mean duration T2. 

Our theoretical approach can be carried over to other 
level configurations and to more than two atoms. For 
given parameters the evaluation should be not too dif- 
ficult. If, however, one is interested in closed algebraic 
expressions the effort will increase considerably with the 
number of atoms. In particular, it would be interesting 
to apply our approach to the situation of the experiment 
of Ref. [0 with its different level configuration and its 
three ions in the trap. 



APPENDIX A: DIPOLE-DIPOLE INTERACTION 
IN THE BLOCK EQUATIONS 

The dipole-dipole interaction enters the Bloch equa- 
tions through r-dependent complex coupling constants 
(cf. Ref. HI) 



^ 2 



1 



ikjir 



(1 - COS^ §j) 



Jk~^~H]^^ (l-Scos^i^,, 



(Al) 



Here "dj is the angle between the transition dipole mo- 
ment Dij and the line connecting the atoms and kji — 
27r/Aji, where Aji is the wavelength of the j-1 transition 
for an atom. For y4.2 ~ one has C2 ~ 0. Thus one 
can neglect the dipole interaction when one atom is in 
state 1 2). The dependence of C3 on r is maximal for 
= 7r/2 and the corresponding C3 is plotted in Fig. 
nfl. For atomic distances greater than about three quar- 




FIG. 10: The complex dipole-dipole coupling constant C3 for 
the strong transition as a function of the atomic distance. 



FIG. 9: The theoretical rate, 712, of double intensity periods 
per unit time shows distance-dependent cooperative effects 
(r23 = 0.5 A3, f22 = 0.01 A3, zero detuning). 

is probably much easier to measure than the double jump 



ters of a wave length of the strong transition, IC3I is less 
than 0.2 A3, but for smaller distances ReC3 approaches 
A3 and ImCs diverges. 

The reset operation TZ and i?cond are given by the same 
expressions as in Rcf. , except for the detuning. One 
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has The summands in Eq. are given by 

7^(p) = (^3 + Rc C3 ) R+ pRl + {A3 - Rc C3 ) R- pRl 

(A2) 

where 

R+ = {3^3 + 3^3) /V2 

= I5>(si3| + |si3)(e3| + (|si2)(s23|-|ai2)(a23|)/V2, 

R- - (5r3-^2i)/V2 

= 1.9) (ai3 1 + |ai3) (63! + (|si2) (023! + |ai2) (S23I) /i^) 
I 



ttO 
cond 



h r 



A3 (|S23)(S23| + |a23)(a23|) + (^3 + C3) |si3) (si3 | + {A3 - C3) |ai3) (ai3 | + 2A3 |e3)(e3| 
V2^]3(|.9)(sl3| + |sl3>(e3|)+^^3(|sl2) (S23I - |ai2)(a23|) + H.c. 



2|e2)(e2| + |si2)(si2| + |ai2)(ai2| + |s23)(s23| + |a23)(a23| 



(A4) 



^^cond(f^2) = ^ V2fl2{\g){.S,2\ + |si2)(e2|) + f^2 (^13) (^23 I + 1 013) (023 | ) + H.C.' 



(A5) 



From Eq. (A4) one sees that ReC3 changes the spon- 
taneous decay rates and that ImCs leads to level shifts. 
Therefore, for small r, the decay rate of |ai3) approaches 
in this case and the large level shifts cause a decrease 
of fluorescence associated with the levels |si3) and |ai3). 



APPENDIX B: CALCULATION OF p{to + At) TO 
FIRST ORDER IN fia 

We write the Bloch equations of Eq. (^) in the form 
p = Cp (Bl) 



where the Liouvillean C = C{A3,il3, A2, C3, ^^2), 
super-operator, can be read off from Eqs. (||) and ( A2) 
(A5). One can decompose £ as 



£ = £0 + >Cr 



where Cq = £(^3, 1^3, A2, C3, 0) and £0^^ = 
-'^{Hcond{^2),p]/h. We note that H^^^^{n2) is Hermi- 
tian and that £0 can be considered as a Liouvillean of 
Bloch equations. Choosing an initial density matrix p(to) 
lying in one of the subspaces in Eqs. (||) - (^) one obtains, 
to first order in 



p{to + At) 



^"^'p{to) 

At 



e 



+ f dTe^°('^*-")£n,e^°" 
Jo 



p{to), (B3) 



just as with usual quantum mechanical perturbation the- 
ory in the interaction picture. Now we use the fact that 
£0, as a Liouvillean of Bloch equations, has an eigenvalue 
(corresponding to steady states) and eigenvalues with 
negative real parts of the order of Q3 and A3. Therefore, 
if At satisfies Eq. (^, the first term on the right-hand 
side of Eq. ( p3| ) gives one of the equilibrium states, 
of £0 given in Eqs. ( |l^ ) - (|l6|), to high accuracy, while 
the term e'~°'^ p{tQ) under the integrand also rapidly ap- 
proaches p^ . After a change of integration variable one 
therefore has to first order in VI2 



p{to + At) ^ p'> + / dTe^°^Cn,p' 



At 



(B4) 



It can be shown that Cn^^p'^ has no components in the 
zero-eigenvalue subspace of £0 Q. Therefore, the in- 
tegrand in Eq. (B4) is rapidly damped, and since At ^ 



(B2) 5I3 , , the upper integration limit can be extended to 



infinity. Hence we can write, to first order in ^2, 



p{tQ + At) = p" + 



POO 

/ dTe^"^£n,P° ■ (B5) 
Jo 



Thus, if At satisfies Eq. (0) then, to first order in ^2, 
p{tQ + At) is independent of At, and one has 

p(to + At) = pO + (e-£o)-i£o.P° (B6) 



to first order in f22 , where the limit e — )■ 
Multiplying this by £ — e gives 

£p(to+At) = £n,(e-£o)"'£n.p" 



-0 is understood. 



o{nl) 



(B7) 
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which is Eq. (|11|). That the transition rates are inde- 
pendent of the particular choice of At follows from Eqs. 
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[60] If£n2P°had a zero-eigenvalue component this would lead 
to a nonunique solution of Eq. (^), which is not the case. 



